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Abstract: The analytic continuation to an imaginary velocity i£ of the canonical partition 
function of a thermal system expressed in a moving frame has a natural implementation in 
the Euclidean path- integral formulation in terms of shifted boundary conditions. Writing 
the Boltzmann factor as exp[— Lq(H — i£ ■ P)], the Poincare invariance underlying a rela- 
tivistic theory implies a dependence of the free-energy on Lq and the shift £ only through 
the combination (3 = Lq\J\ + £ 2 . This in turn implies a set of Ward identities, some of 
which were previously derived by us, among the correlators of the energy-momentum ten- 
sor. In the infinite-volume limit they lead to relations among the cumulants of the total 
energy distribution and those of the momentum, i.e. they connect the energy and the mo- 
mentum distributions in the canonical ensemble. In finite volume the Poincare symmetry 
translates into exact relations among partition functions and correlation functions defined 
with different sets of (generalized) periodic boundary conditions. They have interesting 
applications in lattice field theory. In particular, they offer Ward identities to renormalize 
non-perturbatively the energy-momentum tensor and novel ways to compute thermody- 
namic potentials. At fixed bare parameters they also provide a simple method to vary the 
temperature in much smaller steps than with the standard procedure. 
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1. Introduction 

It is a recurring theme in quantum field theory that a symmetry has far reaching conse- 
quences even when it is softly broken 1 . In this paper we show that, for a relativistic theory 
set up on a space with one or more dimensions of finite length and (generalized) periodic 
boundary conditions 2 , the underlying Lorentz symmetry leads to interesting consequences. 

1 Here with "softly" we refer to any breaking which does not modify the renormalization pattern of the 
theory, e.g. mass terms and (generalized) periodic boundary conditions. 
2 Translational invariance is thus preserved. 
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When continued analytically to an imaginary velocity vector v = i£, where £ 6 M 3 , the 
canonical partition function of a thermal field theory formulated in a moving frame has a 
straightforward definition in the functional integral formalism. It is the ordinary Euclidean 
path integral with shifted boundary conditions in the time-direction [1, 2]. In the zero- 
temperature limit and in presence of a mass gap, the invariance of the theory (and of its 
vacuum) under the Poincare group forces its free-energy to be independent of the shift 

At non-zero temperature the finite time- length Lq breaks the euclideanized Lorentz 
group softly, consequently the free energy depends on the shift (velocity) explicitly but 
only through the combination (3 = Lq y/l + £ 2 . An interesting set of Ward identities (Wis) 
follows. As shown in section 2, they provide a recursion relation among the cumulants of 
the momentum distribution [1, 2], they relate the total energy and momentum distributions 
in the rest-frame, and they suggest new ways to compute thermodynamic potentials. These 
results generalize those found in Ref. [2] to a generic theory and to a generic value of the 
shift For the clarity of the presentation and to avoid unessential technical complications, 
however, we restrict ourselves to bosonic theories in this paper. 

The considerations above extend to a thermal theory set up in a finite spatial volume 
with periodic boundary conditions. The independence of the free energy on the angle 
between the time and the space directions is replaced by relations among partition functions 
of systems with the time and the spatial extensions Lorentz transformed, e.g. 

Z(Lq, Li, L 2 , L 3 ; £i) = Z{ J^ , lJ\ + & L 2 , L 3 ; (1.1) 

with £ = (£i, 0, 0). The Wis are readily extended to finite- volume systems. 

These properties find interesting applications when a theory is discretized on the lattice, 
where a non-zero shift can easily be implemented [1]. Lorentz invariance, which is recovered 
in the continuum limit only, allows one to vary the temperature of the system by changing 
either £ or Lq, i.e. in much smaller steps (at fixed bare parameters) with respect to varying 
Lq alone. Thanks to the misalignment of the lattice axes with respect to the periodic 
directions, the Wis provide new ways to compute thermodynamic potentials numerically 
and new conditions to renormalize non- perturb at ively the energy-momentum tensor. 

2. Thermal field theory in a moving frame 

In this section we focus on properties of a relativistic thermal system in the infinite- volume 
limit. In a moving frame, the total energy and momentum densities are given by ([3], 
paragraph 133) 

e> = -^(e + v 2 p), p' = ^±P- 2 v, (2.1) 
1 — v 1 1 — v 2 

where v is the velocity of the center-of-mass relative to the observer, e and p are the energy 
density and pressure in the rest frame respectively. The enthalpy density (e + p) in the 
rest frame plays the role of the inertial mass density of the system, and its rest volume 
appears contracted by a factor y/l — v 2 in the moving frame. The standard definition of 
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the partition function is 3 ([4], paragraph 2) 

Z(L ,v)^Ti{e- L ^ n - v -^} , (2.2) 

where H and P are the Hamiltonian and the total momentum operator expressed in a 
moving frame. We focus on the Euclidean formulation, where it is natural to continue Z 
to imaginary velocities v = i£ with the Lorentz group replaced by SO(4). The partition 
function 

Z(L , £) = Tr {e-Lotf-it-P)} (2 3) 

corresponds to the ordinary Euclidean path integral with shifted boundary conditions in 
the time-direction [1, 2]. The free-energy density can be defined as usual 

f(L ,Z) = --±-hiZ(L ,£), (2.4) 

LqV 

where V is the volume observed in the moving frame. In the thermodynamic limit the 
invariance of the dynamics under the SO (4) group implies 

/(L ,£) = /(Wl + £ 2 ,0). (2.5) 

In section 4 we derive this equation in the path integral formalism starting from a finite- 
volume system, and provide the functional form for the finite- volume corrections. Eq. (2.5) 
is consistent with modern thermodynamic arguments on the Lorentz transformation of 
the temperature and the free energy [5, 6] (the issue has been debated for a long time, 
see [7] for a recent discussion). Before entering into the details of the derivation it is 
interesting to discuss the origin of this formula, and to anticipate some of the implications 
of the rich kinematics in the boundary conditions that Lorentz symmetry allows for in the 
path-integral formulation. 

2.1 Ward identities for the total energy and momentum 

Relation (2.5) is the source of certain Wis for the energy-momentum tensor, some of which 
were already derived in Ref. [1, 2]. They can be generated in a quasi-automated fashion 
by deriving the free-energy density with respect to L$ and By remembering that the 
cumulants of the total momentum distribution can be written as [1] 

i \rai+n 2 +"3+l rfin x pjln^ x2n 3 

U t , = _/p2"i p2n 2 p2n 3 \ _ V 1 J u u u mm 

«-{2ni,2n 2 ,2n 3 } — y VI r 2 r 3 /c— ^2ni+2n 2 +2n 3 -l Q^m Q^n 2 Q^n 3 J V °' ^ ^_ n ' V 



^=0 



in the thermodynamic limit a plethora of Ward identities among on-shell correlators of the 
total momentum and/or energy are derived by inserting Eq. (2.5) in (2.6). By choosing 
£ = 0, 0}, it is straightforward to derive the master equation 

^-f2n.0.0! , ( 1 d 



v {2nA0} 
-0 



(-ir«(2„-l) !! {i- — }"/(Lo.«)| £rf n = l,2„... (2.7) 



3 We use the notation Lq because this parameter represents the length of the Euclidean time direction 
in the path integral formalism. 
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If we remember that in the Euclidean the momentum operator maps to — > — iTofc, 
where T fMU (xo) = J d 3 iT F (i) with T^ u being the energy- momentum field of the theory, 
an immediate consequence of Eq. (2.7) is the recursion relation (x all different) 

(Toi(xJ) . ..T 01 (x% n )) c = (2n- 1) — { - (Toi(xg) . . .T 01 (x 2 n )) c } n = 2,3. . . (2.8) 

which extends to a generic theory the derivation presented for the scalar one in Ref. [2]. 
Cumulants with non-null indices in the other directions can be related to those in (2.7) by 
cubic symmetry. 

If we define c\ = e — f and recall that the higher cumulants of the total energy 
distribution are given by 

gn-l Qn 



C n ^y(H n )c=(-l) n+1 



n 



f(Lo,£) ^ =Q n = 2,3... , (2.9) 



it is clear that there is a linear relation among c\ , . . . , c n and the n first derivatives of 
the free-energy density. Since Eq. (2.7) gives the k{ 2n ,o,o} as linear combinations of the 
very same derivatives, a linear relation exists among the n first non-trivial cumulants of 
the energy and momentum distributions in the thermodynamic limit. Some details of the 
required combinatorics are summarized in appendix A. The result reads 4 

(2n-l)!! " (2n-i)\ e 
= -j^r ^ (2Lo) c e , (2.10) 

and it shows that the total energy and momentum distributions of a relativistic thermal 
theory are related. Up to n = 4 we obtain 

L o ^{2,0,0} = c i > 

L 3 k {4m = 9ci + 3L c 2 , (2.11) 

Lt k {&m = 225 ci + 90 L c 2 + 15 L\ c 3 , 

L 7 k {8fifi} = 11025 ci + 4725 L c 2 + 1050 L\ c 3 + 105 L§ c 4 . 

As expected c\ = e + p is necessarily positive. Since L\ c 2 is the specific heat, the fourth 
cumulant of the momentum turns out to always be positive. If we remember that in the 
Euclidean (Too) = — e and (Tkk) = P, Eqs. (2.11) can also be written as 

Lq (Toi T i) c = (Too) - (Tn) , 
Ll (T i T i T i T i) c = 9 (T n ) - 9 (T 00 ) + 3 L {T m T m ) c , (2.12) 



where in each correlator the energy-momentum fields are inserted at different times. These 
Wis generalize to all cumulants of a generic field theory those found in Refs. [1, 2]. They 
show that the thermodynamics of a relativistic theory can be studied from its thermal 
momentum distribution and vice-versa. 



'Notice that the coefficients multiplying a are all positive. 
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2.2 Ward identities in presence of a non-zero shift 



When £ 7^ parity is softly broken by the boundary conditions in the compact direction, 
odd derivatives in the do not vanish anymore, and new interesting Wis hold. By deriving 
once with respect to Lq and it is easy to obtain the first non-trivial relation 



6 



l - f 2 



(2.13) 



An interesting consequence of this equation is that the entropy density s of the system at 
the inverse temperature f5 = -Lo\A + £ 2 is given by 



( T 0fc>£ 



(2.14) 



which, by following Refs. [1, 2], can also be written as 

1 d 



lnZ(L ,£) 



(2.15) 



where 7 = 1/a/1 + £ 2 . Ward identities among correlators with more fields can easily be 
obtained by considering higher order derivatives in Lq and For instance by deriving 
two times with respect to the shift components, by using 



Lo(T ok (Lo)0) tc = —(0) $ 

OQk 



(2.16) 



where O is a generic field with support located at a physical distance from the time-slice 
Lq, we obtain 



( T ok)e, 



Lo£k 



^2 ( T 0i T 0j) 



(2.17) 



This equation and Eq. (2.13) can be enforced in regularizations that break translational 
invariance, such as the lattice, to renormalize non-perturbatively the traceless components 
of the energy- momentum tensor, see section 5. By combining Eqs. (2.14) and (2.17), the 
entropy density can also be computed as 



.1 V- ( T 0* T 0j)g,c 

~~ 9 Z_> /T„.\ JT„.\. 



(Toi)^(r 0j -)^ 



and the analogous expression for the specific heat reads 



v ^ _1_ \ ^ 
2 9 Z> 



2 ^ (T 0i )s(T 0j )s e 



2\ x _ n&*i 
V ° £2 



(1 - 2^ - 3 



(2.18) 



(2.19) 
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3. Generalized periodic boundary conditions 

Consider a quantum field theory 5 defined on R d , an orthonormal basis, and d linearly 
independent primitive vectors (/x = 0, 1, . . . , d — 1). The latter can be represented 
by a primitive matrix V G GL(d, M) whose columns are the components of in the 
orthonormal basis. For a given point labeled with the coordinates x^, the field is identified 
at all points with coordinates 

%n + V^m u , m v G Z , (3.1) 

i.e. we impose generalized periodic boundary conditions (GPBCs) 6 . The shifted boundary 
conditions which implement the partition function in Eq. (2.3) are a special case of GPBCs. 

In addition to the parameters already present in infinite volume, the finite-volume 
theory contains d 2 extra parameters specifying the coordinates of the primitive vectors. 
By defining the primitive cell as usual 

Q = [x G R d | x M = V^U, < t„ < l} , (3.2) 

d(d — l)/2 parameters specify the orientation of the cell relative to the orthonormal basis 
(in d = 3 and for orthogonal primitive vectors, these are the Euler angles), while d(d+ 1)/2 
fix its geometry, namely the length of the vectors and the d(d — l)/2 angles between 
them. For a Lorentz-invariant theory, the absolute orientation of the primitive cell is clearly 
of no consequence. The partition function of the finite-volume theory is unchanged if V is 
replaced by 

V -> AV, AG SO(d) , (3.3) 

i.e. the d(d — l)/2 parameters that specify the orientation of the cell are redundant. This 
is the invariance which allows one to generalize Eq. (2.5) in finite volume, and to derive 
the corresponding Wis. At variance with the infinite-volume case, the partition function 
of a finite- volume theory is also left unchanged under the discrete group of transformations 
SL(e£, Z). As is well known from crystallography, two geometrically different primitive cells 
may in fact describe the same crystal: any set of vectors that generates the same 
discrete set of points where the fields are identified is equivalent. This amounts to the 
freedom of replacing the matrix V by a new matrix whose columns are linear combinations 
with integer coefficients of the columns of V, with the restriction that the inverse relation 
exists and contains only integer coefficients. The latter condition requires the determinant 
of the two primitive matrices to be equal up to a sign. Here we restrict ourselves to the 
transformations with positive sign, which maintain the orientation of the unit cell. In short, 
the transformation is 

V — >■ VM, MGSL(d,Z), (3.4) 

i.e. a discrete equivalence between two GPBCs. To summarize, the most general relation 
between two primitive matrices V and W corresponding to a relativistic field theory with 

5 Since the considerations in this section are valid for a generic number of dimensions d > 1, we leave the 
value of d unspecified. In the rest of the paper the results of this section will be used for d = 4. 
6 For brevity we refer to a theory satisfying GPBCs as a finite-volume theory. 
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two different sets of GPBCs and equal partition functions, is given by 



W = A VM, AeSO(d), MeSL(d,Z). (3.5) 

The matrix M modifies the geometry of the primitive cell, while A modifies its orientation. 
The freedom to choose the former is a property of periodic boundary conditions, the free- 
dom to choose the latter is a property of the SO(d) invariance of the infinite- volume field 
theory. The relation (3.5) defines an equivalence relation (in the mathematical sense) be- 
tween the primitive matrices V and W. We will write the relation V ~ W . In Appendix B 
we verify in momentum space that two partition functions defined by path integrals with 
a common Lagrangian density and with equivalent sets of boundary conditions are equal. 
More precisely, the actions in the two theories are related by 

S(V;lcP]) = S(AVM;[<p A ]), (3.6) 

where [4> A ] means that every field of the theory has been rotated. Correlation functions of 
fields can also be mapped between equivalent descriptions of the same system by taking 
into account their transformation properties under the SO(d) group. It is also interesting 
to notice that, by an appropriate field transformation, the effect of the non-orthogonality 
of the original primitive vectors can be re-absorbed into a re-definition of the action, see 
again appendix B. 

4. Finite-volume theory with shifted boundary conditions 

The finite-volume analogue of the partition function 7 (2.3) 

Z(V shc ) = Tr{e' L ^~^}, (4.1) 

where 



Kibe 



/ L 
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Li 













L 2 





\L & 








L ?> ) 



(4.2) 



can be expressed as a Euclidean path integral with the fields satisfying standard periodic 
boundary conditions in the spatial directions, and shifted boundary conditions [1, 8] in 
time 8 

cf>(L ,x) =<f>(0,x-LoZ) . (4.3) 

Due to the spatial periodicity, £' k = + L^/Lq is equivalent to and therefore the 
imaginary velocity components can be restricted to the interval 

-—<&< — ■ (4-4) 
2L " 2L y > 

7 In finite volume we will use the primitive matrix as argument of the partition function Z and of the 
free-energy /. 

8 Relative to these references, we adopt here a different sign convention for the shift in the path integral. 
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For later use it is useful to note the effect of taking derivatives of the partition function 
with respect to the external parameters, 



{Too)v shc 

(Tok)v shc 
i T kk)v A 



1 



L n 







1 d 
LoLiL 2 L 3 d^k 



lnZ(V shc ) 



k = 1,2,3 



(4.5) 



L0L1L2L3 \ "SLfc s,t d£k 
and to introduce the notation 
*2\-V2 



( d d 
L k —+S kl r r ]hiZ(Vs bc ) k = 1,2,3, 



7=(l+* 2 )" 
By defining 



2 , /-2\-l/2 



im = (1 + ^ + ) 



t* = (1 + $ 



2\-l/2 



(4.6) 



Vi = M~ 1 RV shc M 



with 
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( 7i 
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l) 



(4.7) 



(4.8) 



we conclude from section 3 that Z(V S ^, C ) = Z{V\). We first focus on the case £2 = £3 = 0, 
and later use the SO (3) rotation symmetry to generalize the result to a generic shift vector. 
The partition function can be interpreted in terms of the states that propagate in the 
direction given by the first column of V\. In the thermal field theory language, the latter 
are the eigenstates of the 'screening' Hamiltonian H, which acts on states living on a slice of 
dimensions (Lo/71) x L2 x L3 with ordinary periodic boundary conditions. Their spectrum 
yields the spatial correlation lengths of the thermal system at inverse temperature (Lo/71). 
The partition function can thus be written as 



Z(Vi) = Tr { exp -L l7l (# + i^u)}, 



(4.9) 



where u is the momentum operator along the primitive vector of length (Lo/71). Its eigen- 
values are the Matsubara frequencies uj n = 7i^p n£Z. Assuming that the Hamiltonian 
H has a translationally invariant vacuum and a mass gap, the right-hand side of Eq. (4.9) 
becomes insensitive to the phase in the limit Li — > 00 at fixed £1 (with exponentially 
small corrections, see below). This in turn implies that the free energy densities associated 
with V s bc and diag(Li7i, Lo/71, L2, L3) are equal. Thanks to the invariance of the infinite- 
volume theory under three-dimensional rotations, this result extends to a generic imaginary 
velocity £. In the thermodynamic limit the net effect of the generic shift £ is thus to lower 
the temperature from 1/Lq to 1/(3 = 1/(Lq\/1 + £ 2 ), i.e. we have proved Eq. (2.5). As 
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anticipated in section 2, when the primitive cell dimensions are all asymptotically large 
the system is characterized by a single 'short' periodic direction of length (3 = Lo\A + £ 2 
which is interpreted as being its inverse temperature. Its orientation is unusual in that it 
is not aligned along the time-direction, but due to its SO(4) symmetry this is irrelevant. 

In a finite- volume the length of the box dimensions are further sources of SO (4) soft 
breakings, and the above analysis is significantly more involved. To go straight to the 
point, let us assume again that only £i ^ 0. The Euclidean finite- volume counterpart of 
the textbook relations (2.1) read 

{Too)v sbc = 7i ((Too) RV sbc + t,ii T ii)RV shc ) - 2£i7i( 7 oi} J RV sbc , 

((7oo)w sbc - (Tu)Rv sbc + 7i(! - Ci)( T oi)w 6 



( T oi)v sbc 
where in this case 



sbc ' 



(4.10) 



RVsbc 







£i7i£i 





\ 







L171 








?3=0 








L 2 















L 3 J 



(4.11) 



In general the term on the r.h.s (?oi)ijy Bbc does not vanish in finite volume, while it does in 
the thermodynamic limit where Eqs. (2.1) are reproduced. Thus if we consider a thermal 
system satisfying ordinary spatial periodic boundary conditions 'moving' at imaginary 
velocity £, an attempt to 'boost' it back to the rest frame modifies its spatial boundary 
conditions in such a way that the momentum density does not vanish. This effect becomes 
irrelevant when the spatial correlation length is finite and the volume becomes large (see 
below). However, a remarkable property of periodic boundary conditions is that there are 
discrete values of the (imaginary) velocity for which the system at rest still obeys ordinary 
spatial periodic boundary conditions, and the textbook relations in Eqs. (2.1) hold. The 
term (Toi)i?y sbc does vanish in finite volume when the the parameters of the system satisfy 



£i7i 2 6 



= qeZ. (4.12) 
diag(L /7i,Li7i,L 2 ,L 3 ) by an SL(4,Z) 



Indeed Eq. (4.12) implies that i?T4bc|^ 2= ^ 3= o 
transformation, see Eq. (3.4). 

4.1 Finite-size effects 

It is interesting to ask about the magnitude of finite-size corrections to Eq. (2.5). For a 
generic shift the effect of the finite value of L\ in the free energy f{V\) can be quantified 

as 

1 



L0LIL2L3 



£i7i-Evac(^i) - In (l + v e ~ 

l-particle states 



£l7l(£+*£iP 



+ ... 



(4.13) 



where V± is the (0, 0) minor of the matrix V\ which describes the space on which the 
eigenstates of the Hamiltonian H are defined. The vacuum energy -E V ac(^i) on the space 
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V\ corresponds to the free energy of the system in the limit L\ — > oo. As indicated in 
Eq. (4.13), the leading correction 

Xl = v _ y e -ii7i(E+^iPi) (4. 14 ) 

1-particlc states 

comes from one-particle states, where the factor u stands for the multiplicity of the lightest 
screening state of mass M. The allowed momenta in the periodic box described by V± are 
given by 



/2^_ 7i6p2 _ 7i £\ 



p 



L 

P2 



\ P3 J 



p 2 = 2ir T ^, P3 = 2tt T ^, neZ 3 (4.15) 

-^2 -^3 



In the following we assume that a momentum of order 1/Lo always costs a substantial 
gap in energy, and therefore set n\ = 0. Then p is orthogonal to the short periodic 
direction u 1 = (Lo/71, £0^2 > -^0^3)- When the box lengths L2 and L 3 are large, we expect 
the dispersion relation of the one-particle states in momenta orthogonal to u to be the 
ordinary relativistic dispersion relation, due to the emerging SO (3) rotation symmetry in 
the space orthogonal to u. The leading contribution thus reads 

X l = ~ U e ~ Ll 71 ^ +P 2 +Pl +Tl (P2 £.2 +P3?3 ) 2 - i€l 71 (P2 6 +P3 & ) ) _ (4-16) 

L0LIL2L3 ^— * 



P2 ,P3 



Using the Poisson summation formula, diagonalizing the quadratic form under the square 
root and appropriately rescaling the momentum integration variables, we arrive at 

^ V f r -^l7i>/MW+ip-x f417 x 

^ l7 i ^ _y Ha (2tt)2 ^ 

m 2 ,m 3 eZ 

with 



7 (£2 + Cl) L i7igi ~ rn 2 L2^2 ~ m 3 L 3 & m 2 L 2 g3 ~ m 3 L 3 j 2 



(4.18) 



Using the observation that 
72 



f AP e -\ Xl \^M^ + i P - X = _ 2 M dr A 3( M 2 } r = (x 2 + x2) l/2 (41Q) 

7 (27r)^ r 
where A 3 (r, M 2 ) is the propagator of a free massive scalar particle on M 3 , we finally obtain 



1 - 1V V 1 d \ e ~ Mr ] (420) 



where /x = (Li, m2L2,m 3 L 3 ), and Qij = (<5jj + (7 — l)£i£j/£ 2 ) defines a positive norm which 
takes into account the Euclidean version of relativistic length contraction in direction 
The leading contribution to X\ is thus given by the value of (m 2 , 777,3) that minimizes the 
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norm of vector x in Eq. (4.18), i.e. 9 777-2 = mz = 0. The leading finite-volume effect 
associated with L\ is thus 



X = - 7 ^ 
1 27tL L? 7 2 



1 + mLi — 

723 



e _ML l7 / 723 _ ( 4 21 ) 



To obtain the finite volume corrections to E vac (Vi) due to the finiteness of the other spatial 
directions we can proceed iteratively as follows. A helpful observation is that in the limit 
L\ — > 00, the shift — Li7i£i in V\ can be ignored if the ground state is translationally 
invariant. We thus obtain 

7i£vac(F!) = - lim J-lnZ(K; c ), (4.22) 

where V^ hc is obtained from V s b c by making the two-step substitutions: first treat 71 as 
independent of £1 and set 

L Lo/71 , Li L171 , £1 , £ 2 -> 7i6 , 6 -> 7i6 , (4.23) 

and then assign to 71 its value in Eq. (4.6). The partition function in Eq. (4.22) can now 
be interpreted in terms of states living on slices of constant coordinate X2 ■ The associated 
leading finite-size corrections are then given by I2 which, as expected, matches what one 
would obtain by cyclically permuting the three directions in expression (4.21). In summary, 
the leading finite-size contributions to the free energy are 

f(V shc ) - /(Wl + £ 2 ) = Xi +X 2 +X 3 + • • • (4.24) 

and the larger contribution(s) among those on the r.h.s. depends on the particular geometry 
of the shifted boundary conditions. As a test of this result, we perform an independent 
calculation for a free-boson theory in appendix C. It should be noted that, since the leading 
correction arises from one-particle states, the leading finite-size effects are predicted exactly 
by a free-boson theory if one sets its mass to M. The procedure followed in this section 
and the formula (4.24) generalize to shifted boundary conditions those in Ref. [9]. Also in 
this case the leading finite-volume corrections are fully determined once the mass M and 
the multiplicity v of the lightest screening multiplet are known. 

4.2 Ward identities for total energy and momentum 

The equality Z^bc) = Z(Vk), where Vu is defined analogously to V\ for the /c-direction, 
can be used to generate Wis for correlators of the energy and momentum fields in a quasi- 
automated fashion. It suffices to take derivatives with respect to the parameters of the 
primitive vectors. If we derive once with respect to the first WI is given by 

(Tok)v Bbc + ( T ok)v k = , {(Too)v Bbc ~ (Tkk)v sbc } ■ (4.25) 

1 $k 1 Sfc 



9 If L\ 3> £2, £3, several values of m,2 and mz make comparable contributions, and the original repre- 
sentation in terms of a sum over discrete momenta is more useful. However, in that case the finite-volume 
effects associated with the other directions will be the dominant ones anyhow. 
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The second term on the l.h.s proportional to {To k )y k vanishes in the limit L k — > oo, and as 
expected it vanishes also at finite L k if the condition analogous to Eq. (4.12) is satisfied, 
i.e. Lk 2'f k = q £ % ■ By differentiating twice with respect to and by setting ^ = we 
obtain 

Lo (T ok T k))v sbc ,c ~ L k (T ok T ofc ))y sbciC = (T 00 )y sbc - (T kk ) Vshc , (4.26) 
where all insertions in the same correlator are at a physical distance from each other and 




(4.27) 



Analogously the fourth derivative leads to 
Lo {To k T ok T ok T ok )) Vsbc)C - L\ {f ok f ok f ok T ok )) VshctC = 3 {(T 00 )y sbc - {T kk )v shc } + (4.28) 

3 j-kfc {T kk T kk ) VshctC - L (TooT o)y sbc , c | + 6 |Lq {T ok T ok T 00 ) VshctC - L\ (f ok f ok T kk ) VshC!C j 

after some rearrangements of the various terms, having set £ k = again and having inserted 
all fields at a physical distance. This derivation extends to a generic thermal-field theory 
results previously obtained in the scalar field theory, Eqs. (5.3) and (6.15) in Ref. [2]. 
Again, due to the breaking of Lorentz symmetry, Eqs. (4.26) and (4.28) differ from those 
in (2.12) by terms which vanish in the thermodynamic limit. 

5. Applications on the lattice 

The shifted boundary conditions discussed so far provide an interesting formulation to study 
thermal field theories on the lattice. There are many applications that can potentially 
benefit from them. In this section we sketch a few examples with the computation of 
thermodynamic potentials in mind. 

5.1 Renormalization of the energy-momentum tensor 

In the continuum, the charges associated with translational symmetries, i.e. the total en- 
ergy and momentum fields, do not need any ultraviolet renormalization thanks to the Ward 
identities that they satisfy, for a recent discussion see Ref. [2] and references therein. On 
the lattice, however, translational invariance is broken down to a discrete group and the 
standard charge discretizations acquire finite ultraviolet renormalizations. The renormal- 
ization pattern of the energy-momentum tensor depends on the theory under consideration, 
since its field content determines what operators T^ u can mix with. For definiteness the 
discussion below focuses on the SU(iV) Yang-Mills theory, but it applies to the scalar field 
theory as well. 

The energy-momentum field T^ v is a symmetric rank-two tensor. Its traceless part is an 
irreducible representation of the SO (4) group. On the lattice, however, the diagonal and 
off-diagonal components of this multiplet belong to different irreducible representations 
of the hypercubic lattice symmetry group and therefore renormalize in a different way. 
In SU(iV) Yang-Mills theory, they both renormalize multiplicatively. The Wis (4.25) 



- 12 - 



and (4.26) provide two relations among the expectation values of the diagonal and off- 
diagonal components of the energy-momentum tensor. They can be enforced on the lattice 
to compute the overall renormalization constant Zt of the multiplet, and the relative 
normalization z T between the off-diagonal and the diagonal components [10, 11], 



T* = Z T T i, T&-T& = Z T z T (Too-Tu), etc. (5.1) 

where the fields with a superscript 'R' are the renormalized ones. There are many ways to 
implement this strategy in practice. A possible choice is to require a primitive matrix 



( Lq \ 

^ \L n 



~2~ 2^0 

L 



V T = 

\ L J 

for which the condition (4.12) holds, and compute z T as 

_ 3 (Tqi)v t 



(5.2) 



(5.3) 



T 2 (To )v T - {Tu)v T 

while Zt can be determined from (xo / yo, x 2 ^ y 2 ) 

Z T = (Tqo}v t - (T 22 )y T ^ ^ 

z T L (T 02 (x ) T 02 (y)))v T ,c ~ L (T 02 (x 2 ) T 02 (y)))y T , c 

Being fixed by Wis, the finite renormalization constants Zt and z T depend on the bare 
coupling constant only. Up to discretization effects, they are independent of the kinematics 
used to impose them, e.g. the volume, the temperature, the shift parameter, xo etc. Ulti- 
mately which Wis and/or kinematics yield the most accurate results must be investigated 
numerically. 

5.2 Calculation of the entropy and specific heat 

Once the relevant renormalization constants are determined, the entropy density can be 
computed from the expectation value of Tq^ on a lattice with shifted boundary conditions, 

s = z \ T ok)e , 6 r , (5.5) 

by performing simulations at a single inverse temperature value f3 = L$\Jl + £ 2 , and at 
a volume large enough for finite-size effects to be negligible. The latter are exponentially 
small in (ML), where M is the lightest screening mass of the theory. To properly assess 
discretization effects, a set of full-fledged simulations needs to be performed at several 
lattice spacings. A rough idea on their magnitude, however, can be obtained in the non- 
interacting limit of the theory. For the SU(iV) Yang-Mills theory discretized with the 
Wilson action and for the 'clover' form of the lattice field strength tensor [12], discretization 
effects turn out to be rather small, see Fig. 1 for the choice corresponding to £ = (1, 0, 0). 
The details of the calculation are given in appendix D. Once the entropy has been computed 
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Figure 1: Entropy density at finite lattice spacing for the SU(-/V) Yang-Mills theory in the non- 
interacting limit calculated via Eq. (5.5) and normalized to its continuum value ssb/3 3 = 4?r ^ . 
The discretization used is the Wilson action and the 'clover' form of the lattice field strength tensor, 
see appendix D. The inverse temperature is (3 = + £ 2 , and a is the lattice spacing. 



at various values of /3, the pressure can be computed by integrating s in the temperature. 
The ambiguity left due to the integration constant is consistent with the fact that p is 
defined up to an arbitrary additive renormalization constant. 

The entropy density could also be computed directly from Eq. (2.18) without the need 
for fixing the multiplicative renormalization constant. This would require, however, the 
computation of the two-point correlation functions in a large volume. The latter can also 
be used to access the specific heat of the system. From Eqs. (2.18) and (2.19), by choosing 
all Lfc and all equal, for instance, the speed of sound c s is given by 

1 Cy 3 (ToiT 02 )^ c + e(T iT 01 )^ c 

cl s e (ToiT 02 ) $c -(T Q1 T ol )^ c ' 1 • } 

where as usual in each correlator all fields are inserted at physical distance. Note that 
without shifted boundary conditions, the specific heat would require the computation of a 
four-point function of Tofc [1, 2]. Note also that all the computational strategies sketched 
in this section use correlation functions of local operators that require at most an overall 
renormalization constant. The latter can be fixed by Wis in finite volume as described in 
the previous sub-section, and no ultraviolet power-divergent subtractions are needed. 

5.3 The integral method at fixed shift 

Calculations of thermodynamic quantities in lattice gauge theories (see [13] for a recent 
review) usually focus on obtaining the pressure p. In the thermodynamic limit and for a 
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Figure 2: Pressure at finite lattice spacing for the SU(iV) Yang-Mills theory in the non-interacting 
limit. The discretization used is the Wilson action and the 'clover' form of the lattice field strength 
tensor. The inverse temperature is given by f3 — + £ 2 , and a is the lattice spacing. 



homogeneous system, the pressure is equal to minus the free energy, p = — /, and all other 
thermodynamic potentials can in principle be derived from it by taking derivatives with 
respective to the temperature. In the original integral method proposal [14], the pressure p 
is computed by carrying out a line integral of the gradient of the logarithm of the partition 
function with respect to the bare parameters. The integrand is expressed as an expectation 
value of derivatives of the action with respect to the bare parameters a = (ai, . . . , a n ). 
Starting the integration from a point where by convention the free energy vanishes, one 
obtains 

where S is the lattice action and (. . .) a is an expectation value taken at the bare parameter 
set at. The integral is done by keeping fixed Lq/cl, and the temperature is changed by 
varying the lattice spacing. The latter is achieved by varying the bare coupling, but in doing 
so the other couplings (if any) must be adjusted if one wants to remain on a line of constant 
physics. Moreover a subtraction of the vacuum contribution must be made in evaluating 
the path- integral expectation value of , which is usually done at a temperature different 
from the target temperature, or even at zero temperature [14, 15]. Thanks to the integral 
method, many results have been obtained on the lattice for thermal gauge theories [13]. 
Due to the vacuum subtraction and to the integral on the bare parameters at constant 
physics, it remains difficult, however, to reach large temperatures or to apply this method 
to regularizations where the tuning of the bare parameters is technically demanding. 

The integral method can also be applied in the presence of shifted boundary conditions. 
For a given inverse temperature f3 = Lqa/1 + £ 2 , the continuum limit f3/a — > oo of p/3 4 
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Figure 3: Inverse temperature values that become accessible with the use of shifted boundary 
conditions at a fixed lattice spacing a and for different values of L /a. The inverse temperatures 
accessible with a shift in a single direction, £ = (£i,0, 0), are marked by a double circle. 

can be taken at fixed This means that the angles among the lattice axes and the torus 
directions are kept fixed. The discretization effects on pj3 in the non-interacting limit of 
the SU(iV) Yang-Mills theory are displayed in Fig. 2. The plot shows that, in the free 
theory, they can be drastically reduced by using shifted boundary conditions. It remains 
to be seen whether this fact persists in the interacting theory. 

5.4 Temperature scan at fixed lattice spacing 

The possibility of varying the temperature by changing either L^ja or £ allows for a fine 
scan of the temperature axis at fixed lattice spacing. This is illustrated in Fig. 3, where it 
is also compared with the standard procedure of varying Lq/cl only. This fact may turn out 
to be useful in all those cases where the temperature needs to be changed in small steps, 
e.g. study of phase transitions etc. 

The 'T-integral' method [16] is an approach for computing thermodynamic quantities 
related to the integral method which is however formulated directly in the continuum. Here 
the pressure is computed as the integral 



with the integrand computed by Monte Carlo simulations, and the inverse temperature 
(3 is varied by changing Lq/o, while keeping the bare parameters (bare coupling, quark 
masses,. . . ) of the theory fixed. This method has a number of advantages over the method 
based on Eq. (5.7). The subtraction of the vacuum contribution only requires a single zero- 
temperature simulation, and no tuning of bare parameters is required to perform a scan in 





temperature. A significant drawback however is that for a given, realistic {Lq/o, = 8 . . . 20) 
lattice spacing the temperature can only be varied in rather coarse steps. The shifted 
boundary conditions provide a way to almost completely eliminate this drawback. Since the 
integrand is a Lorentz-scalar, its expectation value in the presence of the shifted boundary 
condition is equal, up to discretization effects, to its expectation value in the unshifted 
ensemble at inverse temperature (3 = £o\A + £ 2 . The integrand can thus be scanned in 
much finer steps, and the integral can be computed as 



1 



p(L ) = (l + e) 2 P (LoVl + e) + / dy{\ + y) e(L Jl + £' 2 ) - 3p(L J 1 + £' 2 ) 



o 



(5.9) 



The shifted boundary conditions and the associated Wis also suggest a different implemen- 
tation of the method. Thanks to Eq. (5.5), and by remembering that P 2 -§^p = —s, the 
pressure difference at two temperatures and at a given lattice spacing can be computed for 
instance as 



(5.10) 



where £' k = 1 and £' ± stands for the two components orthogonal to the fc-direction. 



6. Conclusions 

Lorentz invariance implies a great degree of redundancy in defining a relativistic thermal 
theory in the Euclidean path-integral formalism. In the thermodynamic limit, the orienta- 
tion of the compact periodic direction with respect to the coordinate axes can be chosen at 
will and only its length is physically relevant. This redundancy in the description implies 
that the total energy and momentum distributions in the canonical ensemble are related. 

For a finite-size system, the lengths of the box dimensions break this invariance. The 
orientation of the Matsubara cycle relative to the spatial directions does have effects which, 
however, are exponentially suppressed. In the limit of large spatial volume the latter are 
calculable in terms of the mass and multiplicity of the lightest screening state (s) of the 
theory. Being a soft breaking, the correlation functions of the traceless part of the energy- 
momentum tensor still satisfy exact Ward Identities. 

When the theory is regularized in the ultraviolet on a hypercubic lattice, the latter 
singles out a particular reference frame. The overall orientation of the periodic cycles of the 
finite- volume, finite-temperature system with respect to this preferred coordinate system 
affects renormalized observables at the level of lattice artifacts. As the cutoff is removed, 
the artifacts are suppressed by a power of the lattice spacing. 

The shifted boundary condition introduced in [1, 2] constitute a particularly interest- 
ing instance of the generalized boundary conditions described in section 3. In the language 
of the canonical formalism, the energy eigenstates acquire a phase proportional to their 
momentum. This different but equivalent point of view implies that thermodynamic po- 
tentials can be directly inferred from the response of the partition function to the shift in 
the boundary conditions [1, 2], a response which is also encoded in the expectation value 
and in the correlators of the off-diagonal components of the energy-momentum tensor. 
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The flexibility in the lattice formulation added by the introduction of a triplet £ of 
(renormalized) parameters specifying the temporal boundary condition has interesting ap- 
plications. It suggests new and simpler ways to compute thermodynamic potentials, and 
the Ward identities mentioned above can be enforced in a small volume to determine the 
renormalization constants of the energy-momentum tensor components. The temperature 
can be changed either by varying Lq in multiples of the lattice spacing or via the shift 
parameters This results in a much finer scan of its value at fixed bare parameters, a 
feature that may prove particularly useful in investigations of phase transitions (see for 
instance [17]). 
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A. Derivation of Eq. (2.10) 

The Eq. (2.9), which expresses a linear relation between the ci,...,c n and the n first 
derivatives of the free energy density, is readily inverted 

Lo +1 ^ = (-l)" +1 n! £ - L\ c k , n = 1, 2, . . . (A.l) 
k=l 

On the other hand, after some algebra, it is possible to show that 

i ay i ^ (-ir* (2n- fc -i)! h a>f 

LodLo) J - L l n ^ 2-* {k-l)\{n-k)\ °5L§' n " 1 '^--- ^ 

Using first Eq. (A.2) and then (A.l), the derivatives of the free energy in expression (2.7) 
can be replaced by the c n . One then arrives at the desired relation between the cumulants 
of the momentum and the energy operator, Eq. (2.10). 



B. Momentum-space analysis of GPBCs 

Since translational invariance is left unbroken by the boundary conditions, we can expand 
the fields in Fourier modes. The set of momenta compatible with the boundary conditions 
is 10 

]> = [p e R d | Vii = 2TT{V~ l ) v ^n v , n v G z} (B.l) 
10 In crystallographic terminology this is the reciprocal lattice. 
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and the plane wave expansion reads 

faix) = MP) e ip - x = 4>a{x + Vm) , (Vm)„ = V^m u , m v £ Z . (B.2) 
perv 

Clearly, we have the property 

T V = T VM , MeSL(d,Z). (B.3) 

It is instructive to see how the equivalence (3.5) shows up in momentum space. First, 
consider a term in the action of the infinite- volume theory 

r f n jd i \ n 
S?M)=] ( II (2^5 J c n- an ■ ■ ->P n ) <M^) ■ ■■^n(P n ) (2vr) d <5 (d) (E^)' (B ' 4) 

where as usual the delta function enforces momentum conservation in all d directions. Since 
the fields are irreducible representations, Lorentz symmetry is encoded in the property 

C-^ip 1 ,...^)^ 1 )...^) (B.5) 

. — 1„1 A — 1 TL \ 7 /A— 1„1\ I /* - l„n\ 



A-y, . . . , a-vj^(a-V) • • • ^„(A"V), 

i.e. the action density at momenta (pi, . . . ,p n ) of the rotated field 

^(p)^(AW MA _1 p) (B.6) 

is the same as the action density of the original field at momenta {Ar l p\, . . . , hr l p n ). This 
property guarantees in particular that the infinite- volume action of a rotated field is equal 
to the action of the original field. 

In the finite- volume theory, the same form of the action holds, but the integral over 
momenta is replaced by a sum over the set Ty of momenta compatible with the periodicity 
of the field. We can write this contribution to the action as 

Sn(V; [<f>]) = E c n- an (P 1 ' ■■■^P n ) fe^ ^) ■ ■ ■ ^(Pl , (B.7) 

where Vd is the volume of the primitive cell. Clearly Eq. (B.3) implies that the action of 
the two systems parameterized by V and VM are equal for the same field, 

s n (y-M) = s n (VM-[4>})- (b.8) 

Second, we can also write 

Sn{VM)=^i E C" CT "(A-y,..,A-V) 5 Er=ipi ^ 1 (A-V)..0 CTn (A-V), 
d p*er AV 

(B.9) 

and by using Eq. (B.5) one finds 

S n (V;[<j>]) = S n (AV;[<p A ]). (B.10) 

One immediate implication of Eq. (B.8) and (B.10) is that the partition functions of the 
same field theory with equivalent sets of boundary conditions V and W = AVM are equal 

V~W Z(V) = Z(W). (B.ll) 
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B.l An alternative representation of a field theory with GPBCs 

We return briefly to the plane wave expansion of the field, Eq. (B.2), to mention an 
alternative representation of a field theory with GPBCs. The rotation matrix A can always 
be chosen such that V is triangular. By denoting = L^, we can write V = (R^D, 
where R is triangular and all its diagonal elements are unity, and D = diag(Lo, . . . , Ai-i). 
Consider then the field transformation 

fax) = <f>a(R~ U x) • (B.12) 
By expanding 4>{x) in Fourier modes 

fax)= Y, kip)e ip - x , (B.13) 

p=2nD~ 1 n 

in momentum space Eq. (B.12) becomes 

fap) = faRp) . (B.14) 

From both Eq. (B.12) and Eq. (B.13), it is clear that <p fulfills ordinary periodic boundary 
conditions on a primitive cell with d orthogonal sides of lengths (Lq, . . . , L^-i). The effect 
of the non-orthogonality of the original primitive vectors is absorbed into the action for (f>, 

fav; [4>]) = s n (V; [</>}) =^r £ C^CV, • • • ,Rp n ) ^up> faip l )-fa(p n )- 

d p 4 er D 

(B.15) 

The kinetic term of a scalar field theory, for instance, 

c 2 {Rp\Rp 2 ) = -Ip 11 (RiR)p 2 
is a positive-definite quadratic form in this formulation. 



C. Free energy of a non-interacting bosonic theory with shifted boundaries 



In this appendix we compute the free-energy of a non-interacting bosonic field theory in a 
finite volume. This serve to check Eq. (4.24) in the free theory explicitly, and it also shows 
that the latter predicts correctly the leading finite-size effects for a generic thermal theory 
if the mass and the multiplicity are fixed to those of the lightest screening state(s). 

For a generic set of GPBCs, the SO(4) symmetry allows one to cast the primitive 
matrix V in the form 





V 



\z 3 



o\ 



V 



(C.l) 



see Eq. (3.3). The matrix V specifies the spatial periodic directions of the system, and the 
associated three-dimensional reciprocal lattice can be extracted from Eq. (B.l). It is easy 
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to prove that 



/(V)- lim /(V) = ^log(l-e-W^) 

u p 

2 1 e n(-L a u> p +ip-z) 



L det V 9-^0 4-f n 2 



(C.2) 



where w p = \/ p 2 + M 2 . Thanks to the Poisson summation formula 

r =A /n 2 L2 + ( n2:+ vfc) 2 



detV 



where in d-dimensions 



I e ?i(-L(,Wp+ip-zj 



2u>, 



A d (\x\,M 2 ) = J 



ipx 



(2tt)<V + M 2 ' 



(C.3) 



(C.4) 



We can thus write 



/(V) - lim /(V) = E E 

Lo->oo — — 



X |-A 4 (r,M 2 ) 
r or 



Jr=|Vfc| 



= (k ,k) £Z 4 . (C.5) 



By repeating the argument in all fc-directions, i.e. by sending successively Vkk to infinity, 
we arrive at the master equation 



/(V) -/«, = £ 



r <9r 



A 4 (r, M 2 ) 



Jr=|Vfc| 



(C.6) 



where /oo is the free energy of the system on M 4 , i.e. in infinite volume. Since Eq. (C.6) is 
expressed in terms of the norm of all the position vectors equivalent by periodicity to the 
origin, its form is invariant within an equivalence class of primitive matrices. It therefore 
holds for any V G GL(4, R). 

As an application of Eq. (C.6), we consider the case where V is equal to V s h c defined 
in Eq. (4.2), i.e. 



f(V shc )- lim f(V shc ) = J = E 



19 >A, 



r dr 



A 4 (r, M ) 



r=^/n 2 Ll+(nL ^+fj,) 2 



(C.7) 



where fi = (mi-Li, TTI2L2, 7713L3). Expression (C.7) involves 4d propagators, while Eq. (4.20) 
contains 3d propagators. Using again the Poisson formula, one obtains 



/oo 
dne 
-OO 



ld_ 

r dr 



A 4 (r, M 2 ) 



(C.8) 



moSZm^O" 00 

where the argument of the propagator can be rewritten as 

r 2 = ri 2 L 2 + (t/LqS + /x) 2 = ^ (?? L o + 7^ • »f + (Q*i> Qa*) 



(C.9) 
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with Qij being denned below Eq. (4.20). By setting xq = (77L0/7 + 7^ • and by using 
the 'dimensional reduction' relation between the four- and three-dimensional propagators 



<LC[) — A 4 (|x|,M 2 ) 

00 ^0 



dxo 



-^A 3 (r,^ 2 + M 2 ) 
r or 



(CIO) 



- r=|x| 



the following result emerges 



J 



mo m^O 



-i2irmo7 2 (£-fJ,)/Lo 



ld_ 

r dr 



A 3 (r, M 2 + (27rm 0l /L ) 2 ) 



r=\Qn\ 



(C.ll) 



For large spatial box dimensions L^, we can drop all terms but mo = and obtain 



J 



_7_ 

Lr 



id_ 

r dr 



A 3 (r, M 2 ) 



r=\QlA 



+ ... , A 3 (r,M 2 ) = -^-e" 

4irr 



-Mr 



(C.12) 



Taking into account that the relevant terms in fi that give the leading correction in L\ are 
now \i\ = ±Li, /j,2 = H3 = 0, we recover exactly Eq. (4.20) and therefore Eq. (4.24), if we 
set the mass and the multiplicity in the free-theory equal to those of the lightest screening 
state in the interacting theory. 



D. Free bosonic theory on the lattice with shifted boundary conditions 

On a finite- volume lattice specified by the primitive matrix V^bc in Eq. (4.2), the bosonic 
propagator in position space reads 11 

1 L0-1 i(^-p-t)x +ip-x 

= S P 5z4 S m 2 (g-^) + M 2 + 4EL lS m 2 (f)' ^ 

where BZ stands for the Brillouin zone. We are thus interested in the finite sum 

1 L ° _1 e ix (2n£/L -4>) 



and for each value of p we will set 



* = p-€, ^ 2 = M 2 + 4^sin 2 (^) (D.3) 

k=i 

at the end of the calculation. To this end we generalize a well known contour integral 
calculation, see for instance Ref. [18] . The first observation is that 

S (*o) = 7" E 9(e™ e/Lo ^ /2 ,x ) , g(z, x ) = ^ , (D.4) 

Lin ~~ 0J — \Z — Z ) 

u e=o y ' 



lr The lattice spacing is set to a — 1 in this appendix. 
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and using the fact that g(z, xo) = g(—z, xq) we have 

2L„-1 



E(x ) 



M/L -i4>/2 



,x ) 



(D.5) 



The poles of g(z, xq) in the variable z are on the real axis at 



UJ 



- Zl ,... A = ±- ± J (-y + 1 



Consider the integral 



'all 



dz g(z,x ) 



e iL <t> z 2L _ J ' 



(D.6) 



(D.7) 



where the contour r a u contains all the singularities of the integrand. The latter are all 
simple poles located at 



zi, z 2 , z 3 , z 4 ; z e = e 



iire/L -i<t>/2 



* = 0,1,...(2L -1) 



(D.8) 



The integral 7 a n vanishes since the integrand falls off as \z\~ 5 when \z\ — > oo. As a conse- 
quence the sum of all residues is null, and the sum of the residues at Zi equals minus the 
sum of the residues at zg. Since 



z 2L -l +e ine/L z 2L -2 + ^ + e ine(2L -2)/L Z(j) + e in(2L -l)/L 



e iL 4> z 2L _ 1 = ^_ e in£/L 

with z^ = ze 1 ^/ 2 , the residue of the integrand in Eq. (D.7) at Z£ is 



(D.9) 



2L 



i-Kt/L -i4>/2 



Comparing with Eq. (D.5) we have 

dz g{z,x ) 



P z e iLo ^z 2L ° — 1 



,x ) 



1 



dz g(z, x ) 

- z e iL cp z 2L _ i 



(D.10) 



(D.ll) 



where the contour f encircles the poles Z£ but not the poles Zi, while for T it is the other 
way around. Since 



Res 



1 1 

Z UJ 2 — (z — z^ 1 ) 2 



2(1/*? " zf) 



2\ ' 



then 



£(x ) 



1 



E 



-2.T, 



1 



2 ^ elLo4>z 2L _ 1 1/z 2 _ z 

By setting ui = 2sinh(a)/2), it then follows that zf = e ±U) and 

1 



(D.12) 



(D.13) 



S(x ) 



2 sinh cj 



LVXQ 



oiLo4>+Lou) _ I e iL(,4>—L()U> _ 



(D.14) 



- 23 - 



The real and imaginary parts read 

ReEfz ) = sinh ( L o^/ 2 ) cosh[q)(Lo/2 - x )} - sin 2 (L o 0/2) sinh(£x ) 

sinha) (cosh(Louj) — cos(Lq4>)) 

j w \ -sin(Lo^) sinh(wxo) , 

ImSlXo = t — ; — 7~7X7 — 7 1 77 ~x 77 — 777 ■ (D.16) 

2smn(wj (cosn(Low) — cos(LQ(f>)) 

For </> = 0, corresponding to periodic boundary conditions, one recovers the known (real) 
result 

vy_n_ 1 cosh[£(L /2 - gg)] 

MX ° j " 2^mh^ sinh(L ^/2) ' (D ' 17j 

Finally the propagator is obtained by inserting Eq. (D.14) in Eq. (D.l) for each value of p 
after having made the substitutions in Eq. (D.3). 

D.l Expectation value of To^ for the SU(iV) gauge theory 

We are interested in the expectation value of the momentum density operator in the non- 
interacting limit of the SU(iV) gauge theory in presence of shifted boundary conditions. 
We discretize Tqi using the 'clover' discretization of the field strength tensor as described 
in [19, 20]. Using the perturbative expansion and taking the infinite volume limit we 
obtain [20] 

1 , T . l_ v^ 1 f _fp_ sin(27rl/L - pj) cos 2 (;p 2 /2) sin( Pl ) 
2 (^ 2 " IK 01 > L ^ J BZ ( 2vr )3 4gin2( g _ BA ) + 4Z l =iSin 2 (pk/2) 

= 1 77^3 cos 2 (p 2 /2) sin( Pl ) ImS(l) , (D.18) 

Jbz { 27T r 

where in the latter equation for S(l) we use Eq. (D.16) with J 2 = ^Y%=\ srn2 (Pfc/2) and 
(f> = p ■ £. The three-dimensional integral in Eq. (D.18) can be evaluated numerically 
leading to Fig. 1. 
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